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ABSTRACT 

Type lax supernovae (SNe lax) are proposed as one new sub-class of SNe la since they present 
observational properties that are sufficiently distinct from the bulk of SNe la. SNe lax are the 
most common of all types of peculiar SNe by both number and rate, with an estimated rate of 
occurrence of about 5-30% of the total SN la rate. However, the progenitor systems of SNe lax 
are still uncertain. Analyzing pre-explosion images at SN lax positions provides a direct way to 
place strong constraints on the nature of progenitor systems of SNe lax. In this work, we predict 
pre-explosion properties of binary companion stars in a variety of potential progenitor systems by 
performing detailed binary evolution calculations with the one-dimensional stellar evolution code 
STARS. This will be helpful for constraining progenitor systems of SNe lax from their pre-explosion 
observations. With our binary evolution calculations, it is found that the non-degenerate helium (He) 
companion star to both a massive C/0 WD (> 1.1 Mq) and a hybrid C/O/Ne WD can provide 
an explanation for the observations of SN 2012Z-S1, but the hybrid WD-THe star scenario is more 
favorable. 
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1. INTRODUCTION 

Type lax supernovae (SNe lax) were recently proposed 
as one of the largest classes of peculiar SNe (Foley et al. 
2013). Members of this subclass, have several obser¬ 
vational similarities to normal SNe la (Li et al. 2003; 
Branch et al. 2004; Jha et al. 2006; Phillips et al. 2007; 
Foley et al. 2009, 2010, 2013), but also present suffi¬ 
ciently distinct observational properties. SNe lax are 
significantly fainter than normal SNe la (Foley et al. 
2013, 2014). They have a wide range of explosion 
energies (10^®-10®^ erg), ejecta masses (0.15-0.5 Mq), 
and ^®Ni masses (0.003-0.3 M©). The spectra of 
SNe lax are characterized by lower expansion veloci¬ 
ties (2000-8000 km s“^) than those of normal SNe la 
(15,000kms“^) at similar epochs. Moreover, two SNe 
lax (SN 2004cs and SN 2007J) were identified with strong 
He lines in their spectra (Foley et al. 2009, 2013). 

SNe lax are the most common of all types of pecu¬ 
liar SNe by both number and rate, with an estimated 
rate of occurrence of about 5%-30% of the total SN la 
rate (Li et al. 2011; Foley et al. 2013). Most SNe lax are 
observed in late-type, star-forming galaxies (Foley et al. 
2013; Lyman et al. 2013; White et al. 2015). The obser¬ 
vational constraints suggest relatively short delay times 
for the progenitor systems of SNe lax (< 500 Myr, see 
Foley et al. 2014). However, there is at least one lax 
event, SN 2008ge, which was observed in an SO galaxy 
with no signs of star formation (which suggests a long 
delay time, see Foley et al. 2010, 2013). 

It is accepted that SNe lax are either thermonu¬ 
clear explosions of accreting white dwarfs (WDs) or 
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the core-collapse of massive stars. However, the na¬ 
ture of the progenitor systems which give rise to 
SNe lax has not been yet elucidated (Foley et al. 
2013, 2014; Liu et al. 2013a, 2015; McCully et al. 2014; 
Fisher & Jumper 2015). Given the diversity of SNe lax, 
multiple different progenitor channels are likely to con¬ 
tribute to the observed population. The most likely pro¬ 
genitor scenarios currently proposed for SNe lax are sum¬ 
marized as follows. 

1. The Chandrasekhar-mass (Ch-mass) scenario. A 
WD accretes hydrogen (H)- or helium (He)-rich mate¬ 
rial from a non-degenerate companion until its mass ap¬ 
proaches the Ch-mass, Mch ~ 1.4Mq, at which point 
a thermonuclear explosion ensues (Whelan & Iben 1973; 
Nomoto 1982; Nomoto et al. 1984). Here, the companion 
star could be a H-rich donor (either a main-sequence [MS] 
star, a subgiant, or a red giant [RG]) or a He-rich donor 
(He star). It has been suggested that a Ch-mass WD 
can undergo either a deflagration, or a detonation, or a 
delayed detonation to lead to a SN la explosion (Arnett 
1969; Woosley et al. 1986; Khokhlov 1989). In this work, 
we only focus on SNe lax that were recently thought to 
likely be produced from the deflagration explosion Q2 of 
a WD. Turbulent deflagrations in WDs can cause strong 
mixing of the SN ejecta (e.g. Ropke 2005), and the small 
amount of kinetic energy released from deflagration ex¬ 
plosions is in good agreement with the low expansion ve¬ 
locities of SNe lax (Foley et al. 2013). Furthermore, the 
(weak) pure deflagration explosion of a Ch-mass WD has 
recently been proposed as a promising scenario for SNe 
lax (Kromer et al. 2013; Fink et al. 2014; Kromer et al. 
2015). In this specific deflagration explosion sce¬ 
nario, a weak deflagration explosion fails to explode 
the entire WD, leaving behind a polluted WD remnant 
star (Jordan et al. 2012; Kromer et al. 2013; Fink et al. 
2014). Hydrodynamical simulations have shown that 
the off-centerignited weak deflagrations of Ch-mass 
C/0 or hybrid C/O/Ne WDs (Garcia-Berro et al. 1997; 
Chen et al. 2014; Denissenkov et al. 2015) are able to re- 
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produce the characteristic observational features of SNe 
lax quite well (Jordan et al. 2012; Kromer et al. 2013; 
Fink et al. 2014; Kromer et al. 2015), including the faint 
2008ha-like SNe lax (Kromer et al. 2015). However, it is 
impossible to determine the exact ignition condition of a 
Ch-mass WD in this work because the WD is treated as 
a point mass in our binary calculations (see Section 2). 
We therefore assume that all Ch-mass WDs in our bi¬ 
nary calculations will lead to (weak) pure deflagration 
explosions for SNe lax. 

In addition, the pulsational delayed detonation ex¬ 
plosion of a near Ch-mass WD has been proposed 
to be a possible scenario for the SN lax SN 2012Z 
(Stritzinger et al. 2014) which has been suggested to 
come from a progenitor system consisting of a He-star 
donor and a C/0 WD (McCully et al. 2014). However, 
further and more detailed hydrodynamic and radiative- 
transfer simulations for this pulsational delayed detona¬ 
tion explosion model are still needed. 

2. The sub-Ch-mass scenario (e.g. Iben et al. 1987; 
Livne 1990; Woosley & Weaver 1994; Woosley & Kasen 

2011) . The WD accretes He-rich material from a 
companion star (a non-degenerate He star or a He 
WD) to accumulate a He-layer on its surface. If 
an He-shell accumulation reaches a critical value (~ 
0.02-0.2 M 0 ), a single-detonation (i.e., “.la model,” 
see Bildsten et al. 2007; Shen et al. 2010) or double¬ 
detonation is triggered to cause the thermonuclear ex¬ 
plosion of the sub-Ch-mass WD. However, current 
works show that theoretical spectra and light-curves pre¬ 
dicted from the simulations with a thick He shell (C).l- 
0.2 Mq) do not match the observations of SNe la/lax 
(Kromer et al. 2010; Woosley & Kasen 2011; Sim et al. 

2012) . Therefore, we will assume a thin He shell (~ 
0.05 Mq) in our subsequent calculations for the sub- 
Ch-mass scenario. Recent population synthesis cal¬ 
culations for the sub-Ch-mass double-detonation sce¬ 
nario suggest that the rates and delay-time distribution 
from this scenario seem to be consistent with the ob¬ 
servations of SNe lax (Ruiter et al. 2011; Wang et al. 

2013) . However, current explosion models for the double¬ 
detonation scenario (Bildsten et al. 2007; Fink et al. 
2010; Kromer et al. 2010; Sim et al. 2010; Shen et al. 
2010; Woosley & Kasen 2011) show that this kind of ex¬ 
plosion struggles to reproduce the characteristic features 
of SNe lax, for example, the low ejecta velocity of SNe 
lax, low ®®Ni mass, and strong mixing of their explosion 
ejecta. Also, the sub-Ch-mass explosions are not going to 
undergo a weak pure deflagration as in the case for Ch- 
mass WDs. Nevertheless, further investigations are still 
required and numerous complications remain to be solved 
in such a model (Piersanti et al. 2014; Shen & Moore 
2014; Piro 2015). For instance, Piro (2015) suggests that 
including turbulent mixing in He-accreting WDs can lead 
to more than 50% C/0 in the accreted layer at the time 
of ignition, which probably causes the nucleosynthetic 
products to be different, or similar. It is unclear if any ex¬ 
tended explosion model within this scenario can success¬ 
fully reproduce most observational features of the least 
luminous SNe lax. We therefore also include the sub-Ch- 
mass explosion as a potential scenario for producing SNe 
lax in our presented calculations. 

3. The massive-star core-collapse scenario (e.g. 

Umeda & Nomoto 2005; Valenti et al. 2008; 


Moriya et al. 2010). A core-collapse SN is trig¬ 
gered by the gravitational collapse of the Fe core of 
a massive star (with an initial MS mass above about 
10 Mq. that had its H and a significant amount of 
He stripped from its outer layers. In particular, the 
fallback core-collapse explosions of massive stars have 
been proposed to explain the peculiar SN lax SN 2008ha 
because these specific SN explosions could produce a 
low explosion energy and ^®Ni mass to account for the 
properties of SN 2008ha (Moriya et al. 2010). However, 
this core-collapse scenario seems difficult to explain the 
full diversity of SNe lax. Also, a TP-AGB-like source at 
the position of SN 2008ha has been detected by a recent 
Hubble Space Telescope (HST) analysis (Foley et al. 

2014). This source has been further suggested to be a 
candidate of either the bound remnant of the WD, or 
the companion star of the WD, or the companion of the 
massive-star progenitor (although it is very unlikely, see 
Foley et al. 2014). If future observations confirm that 
this source is indeed a bound remnant of the WD, a 
massive star as the progenitor of SN 2008ha would be 
ruled out. Also, pre-explosion imaging for another SN 
lax, SN 2008ge, has ruled out particularly massive stars 
as potential progenitors (Foley et al. 2010, 2015). 

Analyzing pre-explosion images at the SN position pro¬ 
vides a direct way to constrain SN progenitor systems 
(e.g., see Foley et al. 2014, 2015; McCully et al. 2014). 
To date, no progenitors of normal SNe la have yet been 
directly observed, even for detections of relatively nearby 
SNe la, SN 2011fe and SN 2014J. However, the probable 
progenitor system of a SN lax SN 2012Z (i.e., SN 2012Z- 
Sl) has been recently discovered from pre-explosion US'D 
images (McCully et al. 2014). It is further suggested that 
SN 2012Z had a progenitor system consisting of a He star 
donor and a C/0 WD (McCully et al. 2014). 

In the accreting-WD explosion scenario, the WD can 
only be observed directly in our own Milky Way and sev¬ 
eral very nearby galaxies because the WD would be faint. 
Consequently, the companion stars generally play a ma¬ 
jor role in determining the pre-explosion signatures of 
progenitor systems. It may be possible to determine the 
nature of the companion star of SNe lax, and thus their 
progenitors, by analyzing pre-explosion images at the SN 
positions. In this work, we predict pre-explosion signa¬ 
tures of companion stars for different progenitor systems 
that have been recently proposed for SNe lax. This will 
be very helpful for analyzing future pre-explosion obser¬ 
vations of SNe lax. In Section 2, we describe the nu¬ 
merical method we used in this work. The pre-explosion 
signatures of companion stars are presented in Section 3. 
A comparison between our results and the observation of 
SN 2012Z-S1 is presented in Section 4. The conclusions 
are summarized in Section 5. 

2. NUMERICAL METHOD AND MODELS 

We use the Cambridge stellar evolution code STARS, 
which was originally written by Eggleton (1971, 1972, 
1973) and then updated several times (Han et al. 2000). 
The detailed description of the version we used can 
be found in Han et al. (2000) and Han & Podsiadlowski 
(2004). Generally, Roche-lobe overflow (RLOF) is 
treated by the prescription of Han et al. (2000). All our 
stars have metallicity Z = 0.02 and start in circular or¬ 
bits. For H-rich stars (MS, subgiant, and RG), we set the 


3 


ratio of the typical mixing length (1) to the local pressure 
scale height (iJp), a = l/Hp = 2.0, and the convective 
overshooting parameter, 5ov = 0.12, which roughly cor¬ 
responds to an overshooting length of ~ 0.25 Hp. The 
He-rich star is evolved without enhanced mixing, i.e., the 
convective overshooting parameter, dov = 0. 

In our binary evolution calculations, the accreting 
WDs (i.e., C/0 WDs or hybrid C/O/Ne WDs) are 
treated as a point mass, only detailed structures of the 
mass donor stars are solved in the code. We trace de¬ 
tailed mass-transfer of the binary systems until the WDs 
increase their mass to the critical limit (assuming SN lax 
explosions occur). For different companion star models, 
the mass-growth rate of the WD, Mwd, is defined as 
follows. 

For the H-rich donor star channel^ the mass growth 
rate of the WD, Mwd = f?He-MHe = memWtr\, where 
Afne is the mass-growth rate of the He layer under the 
H-shell burning, and |Mtr| is the mass-transfer rate from 
the non-degenerate H-rich companion star. Also, Tyue 
is the mass-accumulation efficiency for He-shell flashes, 
which is shown below in detail (see also Kato & Hachisu 
2004), and ryn is the mass-accumulation efficiency for H 
burning, which is defined as follows (see also Liu et al. 
2015).4 

f -Mcr.n/I-Mtrl, |Mtr| > Mcr.H, 

m={ 1, Mer,H > \Mtr\ > (1) 

I 0, |Mtr| < |Mcr.H, 

where Mcr,H = 5.0 x 10“'^(1.7/A — 1 )(Mwd/M 0 — 
0.4) M0yr“^ is the critical accretion rate for stable H 
burning, X is the H mass fraction, the mass 

of the accreting WD. Here, the optically thick wind 
(Hachisu et al. 1996) is assumed to blow off all unpro¬ 
cessed material if |Mtr| is greater than Mcr.Hi and the 
lost material is assumed to take away the specific or¬ 
bital angular momentum of the accreting WD. The ef¬ 
fect of magnetic braking is also included by adopting the 
description of angular-momentum loss from Sills et al. 
(2000) . If the C/0 WD increase its mass to near the Ch- 
mass limit (ss 1.4M0), we assume the SN lax explosion 
ensues. 

Very recently, considering the uncertainties of the C- 
burning rate, Chen et al. (2014) have suggested that hy¬ 
brid C/O/Ne WDs as massive as 1.3 M0 can be formed 
if the convective undershooting of the C-burning layer 
is large enough, the carbon fraction below the flame is 
largely reduced. These hybrid C/O/Ne WDs in close 
binary systems can eventually accrete H- or He-rich ma¬ 
terial from their companion stars to reach the Ch-mass 
limit and explode as faint SNe (see Garcia-Berro et al. 
1997; Chen et al. 2014; Denissenkov et al. 2015). In this 
work, this hybrid WD scenario has also been included 
in our detailed binary evolution calculations for the Ch- 
mass scenario. Because all WDs are treated as a point 
mass in our calculations, the C/0 WD and the hybrid 

^ Here, the effect of the accretion disk around a WD is not 
considered because whether a disk instability occurs or not dur¬ 
ing the WD-accreting process is unclear (see Hachisu et al. 2010; 
Kato & Hachisu 2012). 


WD channel are distinguished just based on their ini¬ 
tial WD mass (see Section 4). The initial WD masses 
of M(yd = 0.65,0.7,0.8,0.9,1.0,1.1,1.2 and 1.3 M0 are 
considered for the H-accreting Ch-mass scenario (see also 
Liu et al. 2015). With a low initial WD mass of 0.65 M0, 
only two systems (out of about 120 binary systems that 
are calculated) can lead to SN lax explosions, 0.65 M0 
is thus set to be the minimum WD mass for producing 
SNe lax from this channel. A WD mass of 1.3 M0 is the 
maximum initial mass we assumed for a hybrid C/O/Ne 
WD (see Chen et al. 2014; Denissenkov et al. 2015). 

For the He-rich donor star channel, the mass growth 
rate of the WD, Mwd = i7He|M4| (see also Wang et al. 
2010; Liu et al. 2015), where |M/^| is the mass-transfer 
rate from the He-rich companion star, the rjne is the 
mass-accumulation efficiency for He burning, which is 
set to be 

Mcr.He/|M4|, |M/,| > Mcr.He, 

1, Me,,H > IMt'rl > .2) 

f/He, M,t ^ |M4| ^ Miow, ^ ^ 

1 (no burning), |M/J < Miow, 

where Mcr,He = 7.2 x 1 O“®(Mwd/M 0 — 0.6) M0yr“^ is 
the critical accretion rate for stable He burning; M^t is 
the minimum accretion rate for stable He-shell burning 
(Kato & Hachisu 2004; Piersanti et al. 2014); Mow = 
4 X 10“® M0yr“^ is the minimum accretion rate for weak 
He-shell flashes (Woosley et al. 1986); is obtained 
from linearly interpolated from a grid computed by 
Kato & Hachisu (2004). As above mentioned, if the WD 
increases its mass to the near Ch-mass limit (ss I.4M0), 
we assume the SN explosion ensues. Here, the initial 
WD masses of = 0.865,0.9,1.0,1.1,1.2 and 1.3 M0 
are considered for the He-accreting Ch-mass scenario (see 
also Wang et al. 2010; Liu et al. 2015). The lower mass 
limit of 0.865 M0 is corresponding to the minimum WD 
mass for producing SNe lax from this channel, and a 
upper-limit of 1.3 M0 is the maximum initial mass we 
assumed for the hybrid C/O/Ne WD (see Chen et al. 
2014; Denissenkov et al. 2015). 

For low mass-transfer-rate He-accreting case, |M/^| < 
4 X 10“® M0yr“^, a thick layer of helium is believed to 
grow on the surface of the WD. As a result, once an 
He-shell accumulation reaches a critical value (0.05 M0 
is adopted in this work, see Woosley & Kasen 2011), 
we assume a single detonation or double detonation 
is triggered, the sub-Ch-mass WD explodes as an SN 
lax. When the mass-transfer rate is low, |M4| ^ 1 x 
10“® M0yr“^, it is suggested that the flash when the He 
layer ignites is too weak to initiate a carbon detonation, 
which results in only a single He detonation wave prop¬ 
agating outward (Nomoto 1982). We do not specially 
distinguish this low-mass-transfer case in our binary evo¬ 
lution calculations because the condition of He ignition 
is still uncertain (Bildsten et al. 2007; Shen et al. 2010; 
Woosley & Kasen 2011; Shen & Moore 2014). This is 
different from the setup in Wang et al. (2013) who did 
not include the case of |M(j.| ^ 1 x lO“®M0yr“^ be¬ 
cause they only focused on the double-detonation ex¬ 
plosions. In the present work, the initial WD masses 
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Figure 1. Panel (a): H-R diagram of pre-explosion companions of SNe lax in different progenitor scenarios. In this panel, initial masses 
of WDs in the Ch-mass scenario up to a maximum value of 1.2Mq. Panel (b): similar to panel (a), but I.SMq initial mass WDs are 
also included for the H/He-rich accreting Ch-mass scenario. Panel (c),(d) and (e): similar to panel{b), but for different HST/WFC3 band 
magnitude vs. effective temperature. Panel (f): similar to panel (b), but for colour-magnitude diagram (Mf555w-Mf 814W» Mpsssw)- The 
red error bar shows the HST observation of the SN 2012Z-S1 of McCully et al. (2014). The red circle, blue cross, and black square symbols 
corresponding to the MS, RG, and He star donor Ch-mass channel. The green triangle symbols represent the sub-Ch-mass channel. 
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of = 0.8,0.9,1.0 and 1.1 Mq are adopted for the 
He-accreting sub-Ch-mass scenario. Here, 0.8 Mq corre¬ 
sponds to the prediction for the minimum WD mass for 
carbon burning from recent hydrodynamic simulations. 
The detonation of a C/0 WD may not be triggered for 
lower mass (Sim et al. 2012). We also expect that the 
initial C/0 WD masses are below ~ 1.1 Mq as more 
massive WDs are usually formed consisting of O and Ne, 
i.e., ONe WDs (Siess 2006; Doherty et al. 2015, but see 
Ruiter et al. 2013). 

In the He-accreting sub-Ch-mass scenario, the com¬ 
panion star could also be a He WD. However, the long 
delay times of the Hc-WD donor channel (> 800 Myr, 
see Ruiter et al. 2011) are inconsistent with the ob¬ 
servational constraints on the ages of SN lax progeni¬ 
tors (< 500 Myr; Lyman et al. 2013; Foley et al. 2013; 
White et al. 2015), and this channel is generally pro¬ 
posed to produce calcium-rich SNe such as SN 2005E 
(Foley 2015). Therefore, the He-WD donor scenario is 
not addressed in our calculations. 

3. PRE-EXPLOSION COMPANION STARS 

Using the method described in Section 2, we performed 
detailed calculations for a set of binary systems with var¬ 
ious initial WD masses (M/^q), companion masses (M^) 
and orbital periods (P*) in different progenitor scenarios. 
Consequently, the range of companion stars of progenitor 
properties at the moment of SN lax explosions, e.g., the 
bolometric luminosity (L), effective temperature (Teg), 
photospheric radius (P), orbital velocity (t4rb), and sur¬ 
face gravity (g) are directly determined in our binary cal¬ 
culations with the STARS code. In addition, some proper¬ 
ties of binary progenitor systems when the SN explosion 
occurs can also be calculated, e.g., the orbital periods 
of binary systems, the amount of mass lost during the 
optically thick stellar wind phase, mass-transfer rate at 
the moment of SN explosion, and even equatorial rota¬ 
tional velocities of companion stars (if we assume that 
the companion stars co-rotate with their orbits). All 
these pre-explosion signatures of companion stars and 
binary progenitors can be compared with pre-explosion 
observations of SNe lax to help understand the nature of 
the SN lax progenitor (e.g., McCully et al. 2014). 

Furthermore, to facilitate direct comparison with opti¬ 
cal observations, the bolometric luminosity is converted 
to broad band magnitudes by using the same method as 
Pan et al. (2013) under the assumption that the compan¬ 
ion photosphere emits a blackbody radiation spectrum; 


msx 


-2.51ogio 


/ SA(^BA)dA 
/(/°c/A2)SAdA 



( 3 ) 


where Sx is the sensitivity function of a given filter at 
wavelength A, Bx is the Planck function, d is the distance 
of the star, and /° = 3.631 x 10“^° ergcm“^ s“^ Hz“^ is 
the zero-point value in the AB magnitude system. Here, 
the effective temperature (Tea) and photospheric radius 
(R) of a pre-explosion companion star obtained from our 
binary evolution calculations are needed. Specifically for 
this work, the different filters of the HST/WFC3 system 
and their corresponding sensitivity functions are consid¬ 
ered in this study. The absolute magnitudes are calcu¬ 
lated in the AB magnitude system. 


Figure 1 illustrates the Hertzsprung-Russell (H-R) di¬ 
agram in different representations using the HST wave¬ 
bands. As it is shown, companion stars in the H-rich 
donor and He-rich donor channel are independently lo¬ 
cated in the H-R diagram at the moment of SN explosion. 
We can also clearly distinguish the location of companion 
stars between the He-rich donor Ch-mass channel (i.e., 
C/0 WD-l-He star and hybrid WD-l-He star channel) and 
the He-rich donor sub-Ch-mass channel. This indicates 
that we can probably rule out (or confirm) the proposed 
progenitor scenarios for SNe lax with pre-explosion ob¬ 
servations at the SN lax position. But unfortunately, we 
cannot fully distinguish the hybrid C/O/Ne WDs’ com¬ 
panions and C/0 WDs’ companions in the Ch-mass ex¬ 
plosion scenarios because the accreting WDs are treated 
as a point in our detailed binary evolution calculations 
(but see the discussions in the Section 4). 

4. COMPARISON WITH SN 2012Z 

Recently, a luminous blue source (SN 2012Z-S1) was 
detected in pre-explosion HST images that was coin¬ 
cident with the location of SN 2012Z McCully et al. 
(2014), and the non-degenerate He companion star to a 
C/0 WD (originally proposed for SNe lax by Foley et al. 
2013) was suggested to be the likely explanation of the 
observations (McCully et al. 2014). However, they also 
suggested that the possibility of a massive star and an 
accretion disk around the exploding WD cannot still 
be excluded (see Fig. 2 of McCully et al. 2014). Here, 
we compare our model results with the observations of 
SN 2012Z-S1. 

4.1. Comparison with Detailed Binary Calculations 

A comparison between the results of our binary evolu¬ 
tion calculations for different potential progenitor chan¬ 
nels and the observations of SN 2012Z-S1 are displayed in 
Fig. 1. As shown in this figure, the predictions from the 
Ch-mass scenario with a He star donor are consistent 
with the observations, which agrees with the previous 
suggestion for the progenitor system of SN 2012Z (see 
McCully et al. (2014)). 

However, as previously mentioned, it is impossible to 
distinguish the C/0 WDs and the hybrid WDs in our 
binary calculations because the accreting WD is set up 
as a point mass. To further discuss the differences of 
companion stars in the C/0 WD-THe star and hybrid 
WD-THe star channel, we simply assume the range of 
initial WD mass for the C/0 WD-THe star channel of 
0.865-1.2 Mq, the range of initial WD mass for the hy¬ 
brid WD-THe star channel of 1.1-1.3Mq (left column 
of Fig. 2, see also Chen et al. 2014; Denissenkov et al. 
2015). In Section 4.2, our binary population synthesis 
(BPS) calculations for the hybrid WD-THe star channel 
will obtain that most binary systems in this channel have 
an initial WD mass larger than I.IMq, the binary sys¬ 
tems with a lower initial WD mass are quite rare. In 
addition, only few C/0 WDs are formed with a mass 
exceeding I.IMq. This implies our assumptions on the 
initial masses of C/0 WDs and hybrid WDs are appro¬ 
priate. 

The left column in Fig. 2 shows the locations of com¬ 
panion stars at the moment of SN explosions for the 
C/0 WD-THe star and hybrid WD-THe star Ch-mass 
scenario, respectively. It is shown that the observation 
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Figure 2. Left column: similar to panel (b) of Fig.l, but only for the He-star donor Ch-mass scenario. Right column: the distributions 
(in logarithmic scale) of companion stars in the the plane of (logic Teff, logic T), which are obtained from BPS calculations assuming a 
constant star formation rate of 3.5M0yr~^. Here, the C/O WD-I-He star (top -|- middle row) and hybrid C/O/Ne WD-I-He star (bottom 
row) Ch-mass channel are considered, respectively. To better compare with the observation of the SN 2012Z-S1 (the error in red, see 
McCully et al. 2014), we show the results from the C/O WD-I-He star channel by including (middle row) or excluding (top row) the 
calculations with an initial WD mass of 1.2 Mq. 


of SN 2012Z-S1 is only consistent with the predicted com¬ 
panion locations from the hybrid WD-I-He star scenario 
and the C/O WD-I-He star scenario with an initial C/O 
WD mass of 1.2 Mq. Taking the problem of the origin of 
very massive C/O WDs into account, our detailed binary 
evolution calculations seem to disfavor that SN 2012Z-S1 
is a non-degenerate companion star to a C/O WD, it is 
more likely to be a He star with a hybrid C/O/Ne WD. 
Because the hybrid WDs have much lower C to O abun¬ 
dance ratios at the moment of the explosive C ignition 


than their pure C/O counterparts (Denissenkov et al. 
2015), which probably lead to different observational 
characteristics from those of the Ch-mass C/O WDs af¬ 
ter the SN explosions and thus being distinguished by 
spectroscopy observations. Recently, an off-centre defla¬ 
gration in a near Ch-mass hybrid C/O/Ne WD has been 
simulated by Kromer et al. (2015). This showed that de¬ 
flagrations in near Ch-mass hybrid C/O/Ne WDs can ex¬ 
plain the faint SN lax SN 2008ha (Kromer et al. 2015). 
However, only a simple near Ch-mass hybrid C/O/Ne 
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WD model was used, more future works with consider¬ 
ing different initial conditions are needed to explore the 
deflagration explosions of hybrid WDs. 

4.2. Comparison with BPS Calculations 

Furthermore, to obtain the distributions of properties 
of companion stars of C/0 WD (or hybrid C/O/Ne WD) 
+ He-star Ch-mass scenario at the moment of SN explo¬ 
sion, we have performed a detailed Monte Carlo simula¬ 
tion with a rapid population synthesis code (Hurley et al. 
2000, 2002). The code follows the evolution of binaries 
with their properties being recorded at every step. If a 
binary system evolves to a C/0 WD (or hybrid C/O/Ne 
WD) -I- He-star system, and if the system, at the begin¬ 
ning of the RLOF phase, is located in the SN la produc¬ 
tion regions in the plane of (logioH*, M^) for its 
, where P*, and are, respectively, the orbital 

period, the secondary mass, and the WD’s mass of the 
C/0 WD (or hybrid C/O/Ne WD) -|- He-star system at 
the beginning of the RLOF, we assume that an SN la is 
resulted, and the properties of the WD binary at the mo¬ 
ment of SN explosion are obtained by interpolation in the 
three-dimensional grid (M/y^, logipP*) of the close 
WD binaries obtained in our detailed binary evolution 
calculations. We refer to Han & Podsiadlowski (2004) 
and Liu et al. (2015) for a detailed description about the 
BPS method used here. 

The initial mass function of Miller & Scalo (1979) is 
used in this work. We assumed a circular binary or¬ 
bit and set up a constant initial mass ratio distribu¬ 
tion (i.e., n{q') =constant, see Duchene & Kraus 2013). 
We assumed a constant star formation rate (SFR) of 
3.5M0yr“^. The standard energy equations of Webbink 
(1984) were used to calculate the output of the common 
envelope (CE) phase. Similar to our previous studies 
(e.g., Wang et al. 2010; Liu et al. 2015), we use a single 
free parameter oceA to describe the CE ejection process, 
and we set acEA=1.5. Here, oce is the CE ejection effi¬ 
ciency, i.e. the fraction of the released orbital energy used 
to eject the CE; A is a structure parameter that depends 
on the evolutionary stage of the donor star. Specifically 
for the hybrid WDs-|-He star channel, another parameter 
of the uncertainty of the carbon burning rate (CBR, see 
Chen et al. 2014) is introduced and it is set to 0.1 (see 
also Wang et al. 2014). 

The distributions of companion-star temperatures and 
luminosities at the moment of SN lax explosion from 
our BPS calculations for the C/0 WD-|-He star and hy¬ 
brid WD-f-He star scenario are shown in Fig. 2 (right 
column). If we assume the initial maximum C/0 WD 
mass is 1.1 Mq, we cannot obtain companion stars as 
cool as SN 2012Z-S1 in both our detailed binary cal¬ 
culations and BPS calculations. In addition, although 
SN 2012Z-S1 lies in the regions that can potentially be 
reached by our detailed binary simulations for the 1.2 Mq 
C/0 WD-|-He star model (middle row in Fig. 2) and for 
the hybrid WD-(-He star model (bottom row in Fig. 2), 
SN 2012Z-S1 is still much cooler than predictions from 
our BPS calculations. In fact, no star has been obtained 
from our BPS calculations that lies in the region of the 
location of SN 2012Z-S1. We have also studied how the 
different CE efficiencies affect the results. Even if vari¬ 
ous CE efficiencies (i.e., different values of oceA = 0.5- 


1.5) are adopted, we still cannot obtain the stars as 
cool as SN 2012Z-S1. However, we caution that BPS 
calculations still have some uncertainties (see below). 
Also, Fig. 2 shows that our detailed binary calculations 
can significantly contribute to the number of systems in 
the vicinity of SN 2012Z. The discrepancy between the 
STARS results and the BPS predictions is because we do 
not consider the detailed formation channels of WD-|-He 
star binary systems in our binary calculations with the 
STARS code. To obtain a star as cool as SN 2012Z-S1 at 
the moment of SN lax explosion, relatively wide WD-|-He 
star binary systems with an orbital period of > 10 days 
are needed. However, such wide WD-t-He star binary 
systems cannot be produced by our BPS calculations. 
Here, the initial conditions of the populations and the 
treatment of the physics of binary evolution are consid¬ 
ered as reasonable as possible in our BPS studies based 
on current observations and theoretical models (e.g., see 
Duchene & Kraus 2013). If future observations indeed 
confirm that SN 2012Z had a progenitor system which 
contained a helium-star companion, the improvement of 
the BPS method may be required. Taking uncertainties 
on the BPS method into account, we conclude that the 
possibility of a He star donor for SN 2012Z cannot be 
excluded, but the expected probability is very low. 

Additionally, as mentioned by McCully et al. (2014), 
the possibility that SN 2012Z-S1 is a single massive star 
that itself exploded cannot be ruled out yet. Fortunately, 
it may be possible to distinguish between the Ch-mass 
binary model and the single massive-star model through 
HST imaging when SN 2012Z will have faded below 
the brightness of SN 2012Z-S1 (see also McCully et al. 
2014). The simulations for weak deflagration explosions 
of the Ch-mass WD (Jordan et al. 2012; Kromer et al. 
2013; Fink et al. 2014; Kromer et al. 2015) showed that 
an abundance-enriched bound remnant of the WD will 
be left after the SN explosion. Also, the companion 
stars (Marietta et al. 2000; Liu et al. 2012, 2013a,b,c; 
Pan et al. 2013) are expected to survive from the SN 
explosion. If the SN 2012Z produces from the C/0 WD 
(or hybrid C/O/Ne WD)-t-He star scenario, a surviving 
bound remnant of the WD and He companion star prob¬ 
ably be detected in the SN remnants though some dis¬ 
tinct features during their long-term post-explosion evo¬ 
lution such as a spatial velocity (if the binary system 
can be destroyed by the SN explosion, see Jordan et al. 
2012), a particular luminosity evolution (Pan et al. 2013; 
Foley et al. 2014) or even a heavy-element enrichment 
(Liu et al. 2013c; Pan et al. 2013). 

4.3. Uncertainties 

In the SD scenario, only a fairly narrow range in 
the accretion rate will allow stable H- or He-burning 
to be attained on the surface of the WD, avoiding a 
nova explosion. In this work, the prescription of opti¬ 
cally thick wind model from Hachisu et al. (1999) and 
He-retention efficiencies from Kato & Hachisu (2004) are 
used to describe the mass accumulation efficiency of 
accreting WDs. However, strong constraints on the 
uncertainties of the mass-retention efficiencies are still 
lacking (e.g., Shen & Bildsten 2007; Wolf et al. 2013; 
Piersanti et al. 2014). Different mass-retention efficien¬ 
cies are expected to somewhat change the results from 
our binary calculations and thus the BPS results (e.g., 


see Bours et al. 2013; Toonen et al. 2014; Ruiter et al. 
2014; Piersanti et al. 2014). However, we do not expect 
a significant change in the main conclusions presented in 
this paper. This issue will be addressed in detail in a 
forthcoming study. 

In addition, there is a selection effect due to the ini¬ 
tial conditions M2, logioP*) of the populations. 

The initial conditions of BPS calculations may sensitively 
rely on the assumed parameters in specific BPS code, 
which will lead to some uncertainties on the BPS results 
(Toonen et al. 2014). For instance, CE evolution, current 
star-formation rate and initial mass function. However, 
current constraints on these parameters (e.g., the CE ef¬ 
ficiency, see Zorotovic et al. 2010; De Marco et al. 2011; 
Ivanova et al. 2013) are still weak. For a detailed discus¬ 
sion for the effect of different theoretical uncertainties, 
see Claeys et al. (2014). 

5. CONCLUSION AND SUMMARY 

Some potential progenitor scenarios have been pro¬ 
posed to explain the observational properties of SNe lax. 
Although weak deflagrations of Ch-mass WDs seem to be 
the most promising scenario for SNe lax (Jordan et al. 
2012; Kromer et al. 2013; Foley et al. 2013; Fink et al. 
2014; McCully et al. 2014), including the faint SNe lax 
such as 2008ha-like events (Kromer et al. 2015). How¬ 
ever, the nature of the progenitor systems of SNe lax 
is still unclear. In this work, we predict observational 
features of companion stars of binary progenitor systems 
at the moment of the SN lax explosion by performing 
detailed binary evolution calculations with the ID stel¬ 
lar evolution code STARS, which will be very helpful for 
constraining the progenitor systems of SNe lax through 
their pre-explosion imaging (e.g., McCully et al. 2014; 
Foley et al. 2015). 

Comparing our results with the observations of 
SN 2012Z-S1, it is found that our detailed binary evo¬ 
lution calculations for the C/0 WD-|-He star and hy¬ 
brid WD-THe star Ch-mass channel can contribute sig- 
nihcantly to the number of systems in the vicinity of 
SN 2012Z-S1, but it seems the hybrid WD-t-He star Ch- 
mass channel is more favorable for the observations of 
SN 2012Z-S1. However, our BPS calculations for both 
the C/0 WD-THe star and the hybrid WD-fHe star Ch- 
mass channel produces companion stars as cool as the 
SN 2012Z-S1 only at very low probability. We thus con¬ 
clude that the possibility of the He donor star as a com¬ 
panion of SN 2012Z is low based on our BPS predictions, 
but the possibility cannot be excluded. A further con¬ 
firmation needs future observations after SN 2012Z will 
have faded below the brightness of SN 2012Z-S1. 
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